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ABSTRACT 

Galaxy scaling relations, which describe a connection between ostensibly unrelated physical charac- 
teristics of galaxies, testify to an underlying order in galaxy formation that requires understanding. I 
review the development of a scaling relation that 1) unites the well-known Fundamental Plane (FP) 
relation of giant elliptical galaxies and Tully-Fisher (TF) relation of disk galaxies, 2) fits low mass 
spheroidal galaxies, including the ultra-faint satellites of our Galaxy, 3) explains the apparent shift of 
lenticular (SO) galaxies relative to both FP or TF, 3) describes all stellar dynamical systems, including 
systems with no dark matter (stellar clusters), 4) associates explicitly the numerical coefficients that 
account for the apparent "tilt" of the FP away from the direct expectation drawn from the virial 
theorem with systematic variations in the total mass-to-light ratio of galaxies within the half-light 
radius, 5) connects with independent results that demonstrate the robustness of mass estimators when 
applied at the half-light radius, and 6) results in smaller scatter for disk galaxies than the TF relation. 
The relation develops naturally from the virial theorem, but implies the existence of additional galaxy 
formation physics that must now be a focus of galaxy formation studies. More pragmatically, the 
relation provides a lynchpin that can be used to measure distances and galaxy masses. I review two 
applications: 1) the cross-calibration of distance measurement methods, and 2) the determination of 
mass-to-light ratios of simple stellar populations as a function of age, and implications of the latter 
for the stellar initial mass function. 



1. INTRODUCTION 

To develop an understanding of any set of objects, we 
first classify them in the expectation that this will help 
us uncover the rules that describe the set. For stars, this 
systematic approach led to stellar classification, eventu- 
ally to the Hertzprung-Russell diagram, and finally to 
theories of stellar structure and nuclear burning that 
comprise one of astronomy's fundamental successes of 
the previous century. For galaxies, this approach has 
been less successful in uncovering simple intuitive guid- 
ing principles. In part, this failure was due to the ab- 
sence of a comprehensive description of galaxy structure 
akin to that available for stars. Theories of galaxy for- 
mation, currently represe nted mostly by nume rical sim- 
ulations (for example see iSpringel et all I2005T ) . are left 
to describe a loosely tied set of observables that include 
galaxy luminosity functions, clustering properties, color 
distributions, and star formation rates for ensembles of 
systems, rather than the specific characteristics of any 
individual galaxy. As such, even if these models success- 
fully reproduce the existing ensemble observations, our 
understanding of galaxies would be quite different in na- 
ture than our understanding of stars. 

What does this long-running failure to identify sim- 
ple rules of galactic structure signify? Perhaps it re- 
flects a greater underlying complexity to galaxies than 
to stars. Perhaps the formation and evolution of galax- 
ies is so strongly sensitive to different variables that each 
galaxy is an entirely distinct entity and we will never 
find a simple description of galaxy structure that is both 
broadly applicable and sufficiently precise for individual 
galaxies. Seen in this light, it becomes clear that the 
search for the unifying principles of galactic structure is 



in essence an attempt to determine the degree to which 
simple, intuitive guiding principles for galaxy formation 
can accurately describe real galaxies at a non-trivial level 
of detail. At its core, this is an argument between the 
potential value of an analytic description of galaxy for- 
mation vs. the need for numerical simulations. 

In this article, I describe recent work that has demon- 
strated that there is indeed an underlying simple order to 
stellar systems of all types and masses beyond that which 
we can currently explain. The observed low scatter about 
this empirical relationship, which ties together the basic 
measurable properties of galaxies, attests to the presence 
of underlying rules. Perhaps the low scatter arises from 
a galaxy formation version of the central limit theorem 
or, perhaps, it points to a more intuitively meaningful 
connection between the way stars form and are packed 
within dark matter potentials. I frame the discussion of 
the relationship in terms of assumptions and refinements 
to the virial theorem as applied to galaxies to clarify the 
additional constraints and information provided by the 
empirical findings. I will show how this new scaling re- 
lation helps address a number of open questions, includ- 
ing some regarding the nature of SO and low luminosity 
galaxies, while at the same time being applicable to all of 
the galaxies that were well described by the Fundamen- 
tal Plane (FP) and Tully-Fisher (TF) relations. After 
describing this scaling relationship, which is referred to 
as the Fundamental Manifold (FM), in reference to its 
antecedent the FP. I proceed to describe ways in which 
it can be exploited to measure other quantities of core 
astrophysical interest, distances and galaxy masses. 
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2. A LOGICAL PATH TO GALAXY SCALING RELATIONS 

Over time, a disjoint set of rules regarding 
galaxy structure, generally referred to as scaling re- 
lations, have been id e ntified (for som e examples see 
Faber fc Jackie [197^ IK ormendvl 119851: iBurstein et al.l 
19971: iKhosroshahi et al.l 120001: iGrahaml 12002ft . The two 
such relations that are in most common use, and include 
measurements of the internal kine matics of galaxies, are 
the Tully-Fisher relationship (TF; iTullv fc Fishen 1 1977ft 
and t he Fundamental Pla ne (FP; iDiorgovski fc Davisl 
Il987t iDressler et al.l 11987ft . Quantitatively, however, 
these are somewhat arcane parameterizations, with non- 
integer coefficients that are derived empirically and de- 
pend on obse rvational details such as filter passbands (for 
example see iPizagno et al.l 120071 ) . Qualitatively, these 
differ from one another in that they apply to restricted, 
non-overlapping sets of galaxy types and are functions of 
different measured quantities. These relations are there- 
fore at best partial answers in our quest for a compre- 
hensive description of galactic structure. As usual, with 
the benefit of hindsight, one can rework a clearer narra- 
tive. We now trace a straightforward, contiguous path 
to the FP and TF that will explain certain features of 
those relationships, and to the unification of those re- 
lationships into one proposed as being app licable to all 
stellar systems (|Zaritskv et al.l l2006al 120081 ) . 

2.1. Starting From the Virial Theorem 

It is a common misconception that the scaling rela- 
tions are simply rephrased versions of the virial theorem. 
Although their origin lies with that theorem, their exis- 
tence implies additional, non-trivial, physical constraints 
on the nature of galaxy formation. This assertion is clar- 
ified by expressing the virial theorem in a form reminis- 
cent of the FP: 

log r = 2 log V Q - log I - log T + log A - log Bq - C , 

(1) 

where the subscript indicates quantities measured at 
a selected radius, ro: Vq is a measure of the internal 
motions within that radius (typically either the circular 
velocity, velocity dispersion, or some combination), Iq 
is the surface brightness within ro, To is the mass-to- 
light ratio of the matter within ro, Aq and Bq are co- 
efficients arising from the integration of the kinetic and 
potential energy terms in the virial theorem (setting, for 
example, J Q v 2 r 2 dr = AV 2 ), and finally Co is an inte- 
gration constant. In principle, ro should be selected to 
encompass the entire system because the virial theorem 
applies to the system as a whole. However, in practice, 
because galaxies have no well-defined edges and the low 
surface brightnesses of their outskirts make measuring 
these quantities difficult, ro is selected to be a compro- 
mise radius, such as that which encompasses half the 
light of the system, r^. This is the first step away from 
the formal (physically correct) application of the virial 
theorem. 

Additional assumptions and simplifications are neces- 
sary to apply Eq. Q] to real systems. Because A, B, 
C , and T can vary from system to system, and also for 
different enclosed radii within the same system, there 
is no a priori assurance that Equation [T] defines a sim- 
ple, limited distribution of galaxies within the (rQ,lQ,Vo) 
space, aka a scaling relation. In other words, solutions 



of Equation [T] exist for any combination of (ro^OiK)) if 
To, Aq, Bq, and Co are unconstrained, and yet galaxies 
do not populate the entire (ro, Iq, Vo)-space. The more 
confined the distribution of galaxies within this space, 
the more restrictive the constraints on the models. The 
value of a scaling relation is that it quantifies the degree 
to which nature is limiting combinations of these param- 
eters, and, by implication, underlying additional physics 
that we have yet to appreciate that lies beyond the virial 
theorem. If galaxies are found only in limited combina- 
tions of ro, Iq, and Vq, then it must also be true that only 
certain combinations of Aq, Bq, Cq and To are allowed. 
Why? 

2.2. A Key Simplification 

I now take a slight detour in our quest for a comprehen- 
sive scaling relation by considering an important result 
regarding the measurement of galaxy masses. The virial 
theorem, and hence Equation U in a slightly different 
guise, is also the primary pathway to galaxy mass de- 
terminations because it can be used to measure T. The 
diffi culty that underlies all such discussions (for examples 
see. |P agdll952l: IBahcall fc Tremaindll98il iHeisler et al.l 
19851: lErickson et al.l Il987t iLittle fc Tremaind 119871 : 
Zaritsky fc Whitell 19941 ) is evident from our expression of 
the virial theorem in Equation [TJ namely the unknown 
numerical values of A, B, and C, hereafter referred to 
as the virial coefficients. For simple geometries these 
coefficients can be evaluated analytically, for example 
the gravitational potential energy for a uniform density 
sphere is 3GM 2 /5R, leading to B = 3G/5 in this par- 
ticular example. However, in reality these coefficients 
are particularly troublesome because we have no way to 
calculate them without having a full knowledge of the 
gravitational potential and the tracer particle distribu- 
tion function and no guarantee that whatever values we 
adopt remain at least roughly constant from system to 
system. 

In the face of such ignorance, one usually adopts 
the simplest possibility — that these values are the 
same from galaxy to galaxy — and then proceeds to 
use simple modeling or back-of-the-envelope arguments 
to obtain numerical values that are inserted into the 
analogs of Eq. Q] (often referred to as mass estima tors : 
IBahcall fc Tremaind H98ll iZaritskv fc Whitd 11994ft . In 
certain cases, additional data, such as measurements of 
the higher order moments of the line of sight velocity 
distribution, constrain the orbits of t he tracer particles 
(Smit h fc Gra^lTiM IDeionghd[l987t IMerrifield fc^Kentl 
1990), providing independent information on the virial 
cocfficiencts. These approaches typically center on 
the application of the J eans equation (for example see 
IMerrifield fc Kent! Il990|) or Schwa r schild modeling (for 
a few examples se e iVandervoortl Il984t iStatlerl Il987t 
iCretton et al.lH999l ) to help in the interpretation of the 
data, but both require data that is far superior than what 
is typically available. Nevertheless, such studies provide 
critical tests for any less sophisticated method, assuming 
that such methods are not inherently doomed by system 
to system variations in the values of the virial coefficients. 
The key prerequisite to the use of mass estimators is es- 
tablishing that variations in the virial coefficients are not 
t he dominant s o urce o f unce rtainty. 

iWalker et all (|2009l ) and iWolf et al.l (|20Tol ) demon- 
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strated, using a range of dynamical models for spheroidal 
galaxies, that the enclosed mass at the half light radius, 
Mh, as estimated from easily measured observables (size, 
luminosity, and velocity dispersion), is robust to the un- 
known details of the internal kinematics and structure 
of the stellar sy stem (robust to ~ 10% or 0.05 dex, 
iWolf et al.|[20Tot) . Those studies were motivated by the 
desire to measure accurate and precise masses for low 
mass stellar systems as tests of hierarchical structure for- 
mation models and dark matter halo profiles^. However, 
in the language of Equation [TJ those studies show that 
the virial coefficients are materially identical from sys- 
tem to system — if one applies the equation at the half- 
light radius, r/,. Correspondingly, they also show that 
if one attempts to use analogs of Eq. Q] for quantities 
measured at radii other than r^ one will realize more 
scatter in the mass estimates. This work demonstrates 
that the success of the FP, and the extended scaling re- 
lation we discuss here, lies in large part to a fortuitous 
choice of rh as the scaling relatio n's fi ducial radius. The 
value of the IWalker et all (|200l and IWolf et al.l (120101) 
work, in the context of the current discussion, is that it 
codified what had generally been assumed without much 
supporting evidence and also highlighted when and how 
the assumption breaks down. 

Using the results of that work, I calculate a numerical 
value for the co mbined quant it y log Ah — log Bh — C'h 
in Equation [TJ IWalker et all (20091 ) find that M h = 
580r/ 1 er^, where the mass is in solar masses, rh is in pc, 
and o v is the line of sight velocity dispersion^] in km s _1 . 
Rewriting this expression in a form similar to Eq. 1 by 
defining 7/j = Lh/irr^ and converting units, results in 

log r h = 2 log o v - log I h - log T h - 0.73 (2) 

comparing Eqns. Q] and [2] highlights the obvious simi- 
larity and leads us to associate o~ v with V for spheroidal 
galaxies and conclude that log Ah — log Bh — Ch = —0.73. 
The derived value of the combined virial coefficients in 
Equation [2] can be tested by comparing mass estimates 
within rh obtained using the more robust methods (Jeans 
or Sch warschild modeling and/or gravitational lensing 
model (jBolton et al.1 120081 )) to those obtained by apply- 
ing Equation [2] to get Th and then mult iplying by the 
lumin osity Ly t . Following that app roach. IZaritskv et al.l 
(2008) independently (prior to the IWalker et al.l (|2009f ) 
study) found a combined value of the virial coefficients 
of —0.75, in what turned out to be excellent agreement 
with Equation [2j Confirmation that the combined virial 
coefficients are roughly the same from system to system 
comes form the low scatter about Equation 2 of real sys- 
tems with independ ently determined values of Th (see 
IZaritskv et aI1l200cl . 

1 Wolf ct al. (2010) argue further that the true half-mass ra- 
dius, rather than th e projected half-mass radius works best, while 
IWalker ct al. (200Sf) present an analysis in p rojected space. The 
two results are consistent and we opt to use the Walker et al. ( 2009) 
result which involves projecting models rather than deprojecting 
the observations. 

2 We explicitly add the subscript v to a to specify that it rep- 
resents a velocity dispersion. However, in cases where we want 
to highlight some other aspect of a, for example that it is mea- 
sured within the half-light radius, r^, we will drop the subscript 
v and replace it with the subscript h, as in <r^. Nevertheless, a 
always refers to the line-of-sight velocity dispersion if it carries any 
subscript. 



Once the system-to-system stability of the virial coeffi- 
cients is confirmed and accepted, the last remaining un- 
known in Eq. [2] is Th- Observationally, the scatter about 
the Fundamental Manifold (FM) is limited to ~ 0.1 dex, 
and is even smaller for subsamples of galaxies drawn 
from the individual studies that comprise the h eteroge- 
nous dataset in that work ()Zaritskv et al.l [2008h . How- 
ever, even if the scatter were zero, it would still be the 
case mathematically, that any combination of (rh, Oh^Ih) 
would be allowed by Equation [5] if Th is unconstrained^- 
The existence of a scaling relation, where galaxies popu- 
late a very limited region of (rh, <Jh, /^)-space also implies 
that Th is constrained. 

For dark matter free systems, the structure is now en- 
tirely defined if one evaluates Th using simple stellar pop- 
ulation models. Because most of the stellar clusters for 
which the set of (rh, Oh, Ih) exist are old (> 10 Gyr), they 
should have nearly the same value of Th (even if they 
are all of the same age, variations in Th will exist due 
to chemical abundance variations and dynamical evolu- 
tion). With Th set to a constant, Equation [2] describes a 
plane in the log — log ah — log I/j-space. Indeed, Milky 
Way globular clusters not only fall onto a plane, they 
fall o nto a line (jPasauiato fc Bertinl l2008; Zar itskv et al.l 
120111 ). which suggests even further constraints on their 
structure, specifically an underlying relationship between 
two of the three measured parameters that removes an 
additional degree of freedorrQ 

2.3. Onward to the Fundamental Plane 

Giant elliptical galaxies empirically obey the FP, which 
has the form 

\ogr h = (31oga v -j\ogI h + 5, (3) 

where f3, 7, and 8 are numerical coefficients. This rela- 
tionship lacks the troublesome Th that is included in Eq. 

Within the framework of Eq. the validity of the FP 
therefore requires that Th oc <7^~^/^ +7 , for those galax- 
ies that occupy the FP. One particular published fit to 
the FP , r h oc ff i.2±o.o7j-o.82±o.02 from ICappellari et al.1 
;2()(H);. therefore implies that 2 - (3 = 0.8 ± 0.07 and 

1 + 7 = —0.18 ± 0.02, or expressed in another form, that 
Th oc Cy- 8 /^ ' 18 for the giant elliptical galaxies that sat- 
isfy the FP. Although alternative fits of the FP exist 
using different samples of ellipticals that are often ob- 
served in different fil t er pass bands (for e xamples see 
Porgensen et~aT1ll996t fBernardi et all 120031 ). they differ 
in detail rather than in spirit. All of these relations im- 
plicitly require that the mass-to-light ratio have power- 
law dependence on o~ v and Ih- 

The relationship between the total mass-to-light ra- 
tio and the structural properties of galaxies must arise 
from physics that dictates how luminous baryons settle 
into dark matter potential wells. As such, it is crit- 
ical to understand if this relationship holds for more 

3 The only obvious physical constraint on is the lower limit 
defined by the mass-to-light ratio of a purely stellar population. 

4 Qualitatively such constraints are not difficult to imagine. For 
example, a cluster with extremely low surface density might not 
have been able to form simply because the cloud from which it 
would have formed from would have been tidally disrupted. Nev- 
ertheless, the quantitative constraint on formation models provided 
by Eq. [2] could be quite challenging and informative 
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Figure 1. 



T in the B-band for local galaxies from 



van dcr Marcl & van Dokkum (2007) using a Jeans equation anal- 
ysis. Here we take V = cr v and have plotted all 62 of their systems. 
The dotted line represents an estimate of T for a purely stellar pop- 
ulation. The solid line repres ents the expectation bas ed on the cr„ 
dependence inferred from thelCappcllari ct al. (20065) FP. 

than giant elliptical galaxies. Using direct measure- 
ments of from Jeans modeling (Figure [JJ data from 
Ivan der Marel &: van Dokkuml2 007) we see both that the 
power law relationship between and cr v holds for giant 
ellipticals consistent with the slope inferred above, 0.8, 
as it must due to their obeying the FP, and that it breaks 
down for low cr v ellipticals. It is not surprising that the 
power-law description breaks down for low a v because 
extrapolating the power law leads to unphysical values 
of Y/i that are smaller than those of a dark- matter-free 
stellar population (the dotted line in the Figure) . On the 
basis of these data, the FP is manifestly only valid for 
galaxies above a threshold a v ~ 100 km sec . The full 
relationship between T/j and (a v , Ih) must therefore be 
more complex than a power-law. In addition to the flat- 
tening at low a v seen in Figure [TJ there is also evidence 
for a t urnover, or deviations from the FP, at larg e values 
of c v (jZaritskv et alJl2008t iBernardi et alJl2011[ ). 

The FP is not an all-inclusive scaling relation, even 
considering only spheroidal galaxies. The relationship 
between Tft and a v must be of higher complexity that 
presu med in the FP (jZaritskv et al.H2008t iTollerud et"afl 
1 2 1 lh . The FP applies only where a power law descrip- 
tion is an acceptable approximation to this more com- 
plicated relationship. Although the simple power-law 
description of Y/j fails, this failure is not a conceptual 
problem because there was never any physical motiva- 
tion for such straightforward behavior. However, can 
be described as a function of a v even when o~ v is < 100 
km sec -1 if one is will to consider a more complex ana- 
lytic description for the structure of galaxies. 

2.4. Incrementally Adding Rotation 

Even ignoring the galaxies with low o~ v in Fig [TJ 
the FP is only valid for galaxies whose stellar dis- 
tribution is dynamically supported by stellar random 
motions. We must therefore search for a description 
that is more broadly applicable than the FP. In cer- 



Figure 2. Deviation from the F P as a function of v r /a v for 
two ga laxy samples, ellipticals from van dcr Marcl & van Dokkuml 
120071 ) (filled red) and SO's from lBedregal et al.l 12009) (open blue). 
As v T /a v increases the deviations become systematic and negative, 
indicating some dependence on v r is needed in the scaling relation. 
The interpretation of this trend is that the FP is not accounting for 
the dynamical support provided by rotation and that this omission 
causes larger deviations as the rotational support becomes more 
important. 



tain early-type galaxies, lenticulars in particular, but 
also lower l uminosity spheroidals (jDavies et H [1983; 
lBenderlll990f ). rotation provides an important additional 
source of dynamical support. This support is evident 
in the classic diagram comparing the ratio of the ro- 
tational velocity, v r , to the veloci ty dispersion, v r /cr v , 
vs. ellipticity ([Davies et all 119831) . but also becomes 
evident when plotting the deviation from the FP (us- 



ing the parameterizatio n oc a„ 



1.2±0.07 r -0.82±0.02 



from 



ICappellari et al.l (120061)1 vs. v r /a v us i ng da ta for el- 
lipticals (Ivan der Marel fc van Dokkuml 12007) and SO's 
(|Bedregal et al.H2006l ) in Figured 

As galaxies become more rotationally supported, the 
use of only a v to measure their dynamical support 
leads to a systematic deviation from the FP. Therefore, 
a natural extension of the FP invokes a combination 
of a v and v r , referred to here as V, to describe this 
support. The commonly suggested combination is ex- 
pressed as V = y/v% /a + a"l , where a is a parameter 
that is determined by the internal structure of the sys- 
tem. Unfortunately, the existing data are insufficient to 
discriminate between choices of a sugges t ed previously 
dBurstein et al.l Il997t iWeiner et al.l 120061 : iKassin et atl 
I2007t iZaritskv et al.l 120081 ). A stand ard choice is a = 2 
alth ough a = 3 is also a cceptable (Wei ner et al.1 120061 ) 
and IZaritskv et all (|2008f ) fit the data to argue for a = 
2.68. The allowed range in values of a is related to the 
unknown nature of the gravitational potential and the 
tracer particle distribution function. Within an isother- 
mal potential, if the orbits are isotropic, then a = 2, 
while for other potentials and orbital anisotropics the 
value of a ^ 2. Here, we will simply adopt a — 2. Al- 
though the value of a is critical for certain applications, 
such as determining whether subtle variations in the 
value of T/j exists between ellipticals and SO's, it is not 
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so for the discussion here. In general, for pressure sup- 
ported systems V reduces to a v , which is what we used 
in our discussion of giant ellipticals, and for purely rota- 
tionally supported galaxies V reduces to v r /^/a, which 
is what we will use for late type galaxies and then only 
affects the normalization of the late-types relative to the 
early types on the FM. We now proceed to discuss the 
inclusion of late-type galaxies in this contexfl 

It is evident from Figure [2] that some of the scatter 
in the FP, even among ellipticals, comes from neglect- 
ing rotational support and that studies that use the FP 
to search for stellar population differences among early- 
types, for example those searching for differences among 
galaxies in different environments, could face systematic 
errors if the degree of rotational support differs in sys- 
tematic ways. This issues could even trickle down to the 
interpretation of studies tha t do not use the FP, but use 
o~ v as an indication of mass (jZhuet al]l20l3) . 

2.5. Onward to and Beyond the Tully-Fisher 
Relationship 

Unlike the FP, the TF does not involve a radial scale 
and it has resisted improvement (i.e. a reduction in 
scatter) via the addition of other structural parameters 
(|Aaronson fc Mouldlll983[ ). The lack of scale dependence 
is a bit puzzling in that one would expect a galaxy that 
is physically twice as large as another, but with the same 
rotation velocity, to have twice as much mass and hence 
twice the luminosity. As such, the existence of the TF 
implies that the larger galaxy must have a commensu- 
rate lower surface brightness (so that it does not have a 
higher total luminosity) — and that this tradeoff must 
be nearly balanced to avoid the expected scale depen- 
dence. Some hints of a radial dependence are appearing 
with larger datasets and better data. While these stud- 
ies suggest that this balance is not exact ()Courteau et al.l 
2007; Sa intonee fc Sp ekkens 20Tl|), we again see that the 
existence of scaling relations requires that the process 
of galaxy formation results in systematic packing of the 
luminous baryons. For disks, analytic formation mod- 
els have long advocated a connection between the spe- 
cific angular mo mentum of the baryons and surface den- 
sity o r disk size (jFall fc Efstathioulll980l:lDalcanton et al.l 
Il997t ). which goes part of the way to explaining this bal- 
ance. Furthermore, systematic behavior between rota- 
tion curves and da rk matter halos has been identified 
(|Salucci et al.ll2007l ). which demonstrates the existence 
of connections between the formation of the entire galaxy 
and the properties of the inner, observable component. 

Because we seek a universal scaling relation, I avoid 
the TF for now and assert that Eq. [2j inclusive of the 
combined kinematic support from v r and o~ v , should ap- 
ply to late type galaxies. I therefore rewrite the FP as an 
expression for the luminosity, L , for c omparison to TF, 
again using the lCappellari et al.l (|2006[ ) parameterization 
of the FP, 

L = 1.2 log V + 0.181ogl ft + r h + Constant (4) 

5 In principle, one can determine the value of a by requiring 
agreement between early and late type galaxies for a single value 
of the constant in Eq. [2] However, in practice because of the 
dichotomy between FP and TF studies, the samples of early and 
late type galaxies have been observed and analyzed differently and 
small differences in photometric systems or analysis technique are 
sufficient to explain the offsets found. 
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Figure 3. Predictions of galaxy luminosity from TF and our re- 
written FP for disk galaxies. We compare the predicted luminosi- 
ties (x-axis) to measured luminosities using the galaxies H ubble 
distance (y-axis) for galaxies from the Pizagno et al. ( 2007) study 
using their best-fit TF r elationship (left panel) to and that using 
our re-expression of the ICap pcllari et al. (2006) FP relation and 
the use of the V parameter (right panel). The scatter is evidently 
smaller in the right panel (0.75 in the left one, 0.47 in the right), 
showing that including a radial scale, as is done in the FM, adds 
significant information. 

In Figure [3l I compare the predictions of L from the TF 
(i.e. using the circul ar velocity and the best fit TF re- 
lation presented by iPizagno et al.l 12007) and from the 
rewritten FP (Equation 2J) with the luminosities inferred 
using distances calculated from their recessional veloci- 
ties for disk galaxies with V > 100 km sec - 1 (because the 
FP ap plies only to such galaxies) . I use the lPizagno et al.l 
(|2007f ) study because it is one of the few TF samples that 
includes measurements of . The scatter is visibly lower 
in the re-written FP relation than in the TF relation 
(0.47 vs. 0.75), demonstrating that the FM parameter- 
ization is at least as good as the best-fit TF. Previous 
studies that searched for a scale dependence may have 
had a difficult time finding one because they either used 
different radii than rh , such as the disk scale length or an 
isophotal radius, or faced additional noise from not using 
the rotational velocity at or near rh- Although for disk 
galaxies, with their mostly circular orbits, one can argue 
that M oc rvf. should be a fairly good approximation re- 
gardless of which of the standard fidu cial radii (disk scale, 
isoph otal, half light) is chosen (see lYegorova fc Saluccl 
120071 for empirical support for this claim), the only one 
of these radii that has a direct connection to the total 
luminosity is the half light radius. 

To see how the TF is a subset of the relationships al- 
lowed by Eq. [2 I begin with the TF relation 



L = Av 5 rl 



(5) 



where A and 5 are constants, rewrite the expression using 
the definition of Ih , and rearrange it to be in the form of 

Eq. m 

\ogrh = 8\ogv r — log Ih + Constant. (6) 

Again there is a strong similarity to Eq. [2] and the 
validity of both equations for giant spiral galaxies sug- 
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rals. Because empirically 6 

l2007tlCourteau et alJ2 007). decreases as v r increases. 
Therefore, it is wrong to accept that the TF applies in 
general because, like the FP, it too faces a fundamental 
conceptual flaw if extended beyond the class of galaxy 
from which it was derived. The FP predicts unphysi- 
cal values of for low <r v systems. The TF predicts 
unphysical values of Y/j for high v r galaxies. Therefore, 
just as in the case of the FP, some curvature in the rela- 
tion is necessary for disk galaxies. Suggestions of nonlin- 
earity in the TF relation have existed for over 20 years 
(jMould et al.lll989t IPersic fc Saluccilll991l) . Neither the 
FP nor TF hold across all galaxy masses, and a more 
complex relation is needed between and V. 

2.6. Closure on a Universal Galaxy Kinematic Scaling 

Relation 

Returning to Equation [TJ we see that the three ingre- 
dients involved in arriving at a universal galaxy scaling 
relation are 1) the virial theorem, 2) the robustness of 
mass estimates at the half light radius, which affirms 
that the virial coefficients are relatively independent of 
the detailed structure and kinematics of each galaxy, and 
3) the relationship between and the directly mea- 
surable structure parameters (rh,Vh, Ih)- The first is a 
natural expectation for collapsed systems, particularly 
sin ce ru <C r V j r i „i . Th e sec ond has bee n dem onstrated 
by IWalker et all IpOOfll ) and IWolf et al.l (poTnh but has 
been evident in the low empirical scatter of the scaling 
relations for decades. The last remains perhaps the most 
mysterious because it involves a connection between the 
efficiency of star formation and the packing of those stars 
within dark matter halos of galaxies. Understanding this 
behavior is particularly challenging given that the val- 
ues of T/j for galaxies vary between values of ~ 1 and 
~ 1000. This behavior lies at the core of understanding 
galaxy formation. 

The classic scaling relations, FP and TF, implicitly 
adopt a power-law scaling between and the observ- 
ables. This description fails as we proceed to either low 
(Figure Q]) or high o~ v galaxies. The power-law assump- 
tion must be abandoned. A more complicated manifold 
relates the observ ables to and was na med the Funda- 
mental Manifold (|Zaritskv et al.l l2006al ) in reference to 
its antecedent the FP. In the interest of completeness, 
I pre sent one speci fi c par ameterization of that surface 
from lZaritskv et al.l (| 20 111 ) 

T h = 1.49 - 0.32 log V h - 0.83 log I h + 0.24 log 2 V h + 

O.Wlog 2 !*- 0.02 log V & log J h . (7) 

However, because there is yet no set of homogeneous data 
that includes internal kinematics and spans all galaxy 
types and luminosities, the derivation of this surface re- 
mains preliminary and varies quantitatively depending 
on what sample is used. Additionally, there is no phys- 
ical motivation for why the surface has any particular 
functional form. As such, larger samples that cover the 
full parameter space are critically needed. Should sub- 
sequent investigations define a completely different func- 
tional form for T/j, that finding would not invalidate any 
of the discussion presented so far other than Equation [7l 
Similar curvature is seen in the behavior of T in 



efforts to match the halo and stellar mass functions 
(|van den Bosch et al.ll200l lYang et al.l I2005T) and lens- 
ing mass measurements and lu minosities ( Hockstr a et al.l 
120051 iMadclba um et al.lfeOOSf ). Because the relationship 
is currently only empirical, the finding that all stellar sys- 
tems can be modeled with a simple, smooth functional 
form is perhaps surprising — but it hints at rules gov- 
erning galaxy formation that are currently not fully un- 
derstood. 

2.7. What Comes Next? 

On the observational front, the refinement of the FM 
requires a homogeneous dataset that includes all galaxy 
types and luminosities and which is volume representa- 
tive (not necessarily complete). Such data would enable 
1) a determination of a, 2) tests of whether the scatter 
increases for certain classes of galaxies, which would ad- 
dress whether the virial coeffciients are sufficiently con- 
stant across galaxy types, and 3) a determination of how 
the surface relating Th to (Vh, Ih) is populated to identify 
additional constraints on galactic structure. Care should 
be taken to minimize sources of photometric scatter by 
choosing passbands that minimize variations arising from 
stellar population variations and extinction. The sample 
should consist of galaxies with independent distance es- 
timators so that distance uncertainties contribute mini- 
mally to rh- Currently, the distances are calculated us- 
ing the recessional velocities and an adopted value of the 
Hubble constant. 2-D kinematics would help address is- 
sues related to inclination corrections for rotation speeds 
and help in the measurement of o~ v and v r in systems 
where both contribute noticeably to the dynamical sup- 
port. The sample might also include more systems where 
Th can be estimated in an independent manner, for ex- 
ample through the use of gravitational lensing. Finally, 
and particularly important for questions related to the 
nature of the luminous baryon distribution, it is neces- 
sary to have independent estimates the stellar mass-to- 
light ratio, T*. One such approach at m easuring T„, 
is to use infrared lumin osities and colors (Mei dt et al.l 
120121 lEskew et al.l 120121 ) , although these estimates de- 
pend on stellar population models, that have their own 
set of uncertainties. Should these uncertainties be sorted 
out in subsequent work, we would be able to use the 
FM and these estimates to uncover the dark matter frac- 
tions for all galaxies within in a way now done only 
with mor e sophisticated modeli ng on smaller samples of 
galaxies (jCappehari et al.ll2012l ). 

3. DISCUSSION 

3.1. The Bifurcation of Stellar Systems 

The function that describes the relation between Th 
and (Vhjh) is roughly parabolic along the Vh axis and 
power-law like along the Ih axis (Equation \7§. The 
general shape places constraints on how baryons set- 
tle and become stars within dark matter potentials. 
However, an important aspect that I have mostly ig- 
nored so far is how this surface is populated. In 
particular, there is a bifurcation in the population 
at low Vh where two branches develop, one head- 
ing upward to large values of T^, populated by the 
Local Group dwarf spheroidals (iZaritskv et all [2006a; 
iTollerud et all 120111: iSalucci et all 120121) and~ultrafaint 
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galaxies (|Zaritskv et al.ll201lHWillman fc Straderll2012l ). 
and another heading downward to low values of T^, pop- 
ulated by ultracompact dwarf galaxies a nd star clusters 
([Mieske et al.ll2008t IZaritskv et all 120111 ). The question 
is whether this bifu rcation points to a problem with 
the scaling relation (jForbes et al.l [2008) or whether it 
is pointing to a fundamental difference between sys- 
tems with and without dark matter halos (Zaritsky et al. 

unm). 

The bifurcation may lie at the heart of a phys- 
ical distinction between galaxies and star clusters 
([Willman fc Straderl [2012] ) and the FM may be a way 
to help explore that division. In particular, it would be 
of interest to determine whether there are objects that lie 
between t hese two branches, or w hether the branches are 
absolute (Misg eld fc Hilkerll2011[ ). Unfortunately, strong 
selection effects come into play. At the extremely low sur- 
face brightnesses of the ultrafaint galaxies, these systems 
require star counts to detect and are therefore not within 
reach if they lie outside the Local Group. It may be pos- 
sible, with the next generation of sky surveys, to detect 
large numbers of such objects eve n though we are con - 
fined to finding only the local ones (jTollerud et al.l f2008). 
On the high surface brightness side, these objects become 
confused with stars dHilker et al.lll999tlDrinkwater et al.l 
l2000t IPhillipps et all 120011) and therefore require unbi- 
ased redshift surveys to identify (for recent examples see 
iMisgeld et aTll201lt IChiboucas et al.ll2011| ). There are a 
handful o f Galactic objects th at might lie between the 
branches (Zari tsky et al.ll2011[ ) but these systems are all 
of low mass, highly susceptible to dynamical effects, and 
may therefore not satisfy the basic requirement of Eq. [TJ 
which is that the system satisfy the virial theorem. 

3.2. Theoretical Work 

The physics of galaxy formation determines the rela- 
tionship between and the structural parameters of 
each galaxy. As such, we desire a theory that will ex- 
plain the principle observables of each galaxy individ- 
ually. Many attempts to reconcile the theory of galaxy 
formation as currently understood to observed scaling re- 
lations (eithe r one for spheroids or disks) exis t (for some 
examples see iDutton fc van den Boschl 120121 and refer- 
ences therein). The key to these models is how one fol- 
lows the complicated physics of galactic feedback and an- 
gular momentum transfer between the baryons and dark 
matter. The models can reproduce general galaxy prop- 
erties if one artificially imposes how these effects scale 
with halo mass. Of c ourse, the next question i s why 
such recaling exists. IDutton fc van den Boschl ([20121 ) 
conclude that certain processes, such as the angular mo- 
mentum evolution of galaxy disks, are not yet sufficiently 
well understood to be modeled sufficiently accurately to 
reproduce the scaling relations. Empirical scaling rela- 
tion provide challenging benchmarks for the models. Ex- 
tending the range of these scaling relations, for example 
down to ultrafaint galaxies, is invaluable because it places 
even stronger constraint s on hypothesized physical m ech- 
anisms. For exa mple, lAnderson &: Bregmanl ([20101 ) and 
iMcGaughl (|2012l ) used the baryonic TF relation to argue 
that there ar e basic concept ual problems with feedback 
models, while IDutton! ([2012D then used this argument to 
examine how to obtain plausible models for star forma- 
tion and feedback. This article is not intended as review 



of theoretical models, but such modeling is critically de- 
pendent on the best possible empirical constraints. Be- 
cause current simulations lack the ability to treat the 
physical processes in detail and realistically, they are not 
reliably predictive and must constantly be compared to 
the available constraints. 

4. APPLICATIONS OF THE FM 

Although understanding galaxy formation remains the 
principal goal behind refining and understanding the FM, 
the relationship also has various uses that do not require 
a theoretical underpinning. For example, as with the 
FP and TF before it, the FM can be used as a distance 
estimator because is in physical (distance dependent) 
units. 

4.1. Using the Scaling Relationships for Distance 
Measurements 

The advantages that the FM provides over its an- 
tecedents come from its universality. For example, the 
FM is applicable to low luminosity local systems, in 
which individual stars are resolved but neither the FP or 
TF apply. This capability is critical because certain dis- 
tance methods require resolved stellar populations, which 
are beyond our current technology's reach for most of the 
galaxies that we know satisfy the FP and TF. As such, 
the FM can be used to cross-calibrate distance estimators 
even when those estimators are not found in the same 
gal axy. For example, Fig ure 4 (reproduced from Figure 
2 of lZaritskv et al.1 ([20111 )) illustrates why it is difficult to 
compare distances obtained from surface brightness fluc- 
tuations (SBF) and Cepheids if one requires having dis- 
tance estimates for the same galaxy from the two meth- 
ods. Likewise, certain methods tend to work over limited 
distance ranges and are therefore not suitable for com- 
parison with other. Using the FM as a fiducial, distances 
obtained using SN la, Cepheids, SBF, the luminosity of 
the tip of the red giant branch (TRGB), circumnuclear 
masers, eclipsing binaries, RR Lyrae stars, and the plan- 
et a^neJoufaeJjuuin^silT functions (PNe) were compared 
bv IZaritskv eTall (f201lh . 

This comparison is done by placing galaxies with in- 
dependent distance measures on the FM, using the dis- 
tance from each particular distance estimator to convert 
rh from angular to physical units. I n Figu re [5] we repro- 
duce Figure 7 from IZaritskv et aD ([20111 ) that presents 
the FM relationship derived using the distance measure- 
ments obtained from each type of distance estimator for 
which there are at least 10 galaxies with all the neces- 
sary data. The velocity dispersion dominated systems, 
which have measured velocity dispersions and sometimes 
also have stellar rotation values, are plotted in the up- 
per panels of the Figure and the purely rotationally sup- 
ported galaxies, which have no quoted stellar velocity 
dispersion and where the inclination-corrected H I width 
is adopted as a measure of v r , are plotted in the lower 
panels. In general, if there is a stellar velocity disper- 
sion measurement, the galaxy is a pressure supported 
systems [y r ja v < 1). Within each panel a color-code 
describes morphology, dividing the early and late type 
galaxies at a T-Type of 1. The zero point of the FM is set 
using the results from the pressure-supported, SBF sam- 
ple, which is the largest subsample of galaxies and also 
predominantly consists of early-type galaxies, which are 
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Figure 4. Distribution of morphological types as a function of 
distance measuring technique among a sample drawn from NED- 
ID. The Y-axis is arbitrarily normalized to match scales among 
subsamples. Differences among the types of galaxies accessible 
to various methods are clear (adopted directly from Figure 2 of 
IZaritskv et~aT]|201ll) 

less susceptible to extinction and stellar population vari- 
ations. Although many more galaxies than those shown 
have distances measured using at least one of these meth- 
ods, many of those lack the measurements necessary to 
place them on the FM. The solid lines represent the FM 
and the only "free" parameter involved here is the nor- 
malization applied to get the mean SBF relation exactly 
(on average) on the FM. Normalization or slope errors for 
any particular distance estimator suggests that there is a 
problem with that estimator. Increased scatter suggests 
a lower precision for that particular distance estimator. 

A cursory examination of the panels reveals no serious 
problems with any of the distance estimators as applied 
to any of the galaxy subsamples, even though some per- 
form better than others (either in terms of zero point or 
scatter). There are a few individual galaxy outliers, al- 
though it is often the same galaxy that is an outlier in 
multiple panels because distances are sometimes avail- 
able from multiple estimators, suggesting that fault lies 
not with the distance estimate but rather with one of the 
other parameters that enters the FM. There are no sta- 
tistically significant differences in zero poi nt between the 
distance estimators (Zarit sky et al.l [201 If) . although the 
allowed differences are still above the level of precision (~ 
percent) that is the goal of current studies. Increasing 
the sample significantly, which is well within what can 
be done with a reasonable investment of resources, will 
provide stricter limits on potential systematic differences 
a mong estimator s. 

Zari tsky et al.l ()201lD also examined the FM residuals 
within a given distance estimator vs. other potential 
sources of systematic error. For example, they examined 
the relationship between FM residual and host galaxy 
metallicity for distances derived using SNe la, they found 
a correlation, confirming previous clai ms of metallicity- 
dependent correction to the la distances dG allaghe r~et al.l 
[2008t iHowell et al.ll200ft iRomaniello et al.ll2008ft ~ 

Once established as a distance estimator itself, the 
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Figure 5. FM using measured distances. A comparison of the 
FM's obtained using different distance estimators. The sample is 
divided into pressure supported (upper panels) and rotationally 
supported (lower panels). Furthermore, color and shape codes dis- 
tinguish galaxy morphologies (blue squares for late type, red cir- 
cles for early type). The x axis is log in kpc, and the y-axis 
is 21og V — log//j — logTj, — C, where V is in km sec -1 , //, is in 
L0/sq. pc, and Tj, is in solar units. C is obtained by calibrating to 
the sample of surface brightness fluctuation distances (SBF) for the 
pressure supported galaxies. The line is the 1:1 e xpectation, not 
a fit to the data (adopted directly from Figure 7 of lZaritskv et all 
120111) . 

FM can be used for a large set of galaxies . In F igure 
I reproduce Figure 12 of IZaritskv et al.l (|2011[ ) that 
shows the relationship between recessional velocity and 
distance obtained from the FM for the same set of galax- 
ies shown in Figured] Comparing between different nor- 
malization of the distances, using the various di stance 
estimators available to them, IZaritskv et al.l ()2011f ) cite a 
systematic uncertainty of 4 km sec -1 Mpc -1 in Hq. The 
inset shows the result of fitting the relationship for galax- 
ies with v > 1500 km sec -1 , 3er outliers excluded, binned 
by 10, with the fit forced through the origin. The low v 
region is excluded to minimize the effect of local flows. 
They find a best fit slope corresponding to Hq = 78 ± 2 
(random) ± 4 (systematic) km sec -1 Mpc -1 . The es- 
timate of the systematic uncertainty does not include a 
variety of potential problems that they ignored (modeling 
of bulk flows, internal extinction corrections, adjustment 
for potential biases in the galaxy sample, etc.), but this 
calculation was done primarily as a plausibility exercise 
to demonstrate the use of the FM. 

One of those sources of uncertainty, internal extinction, 
would be mitiga ted by going to the infrared , as shown 
most recently bv lFreedman fc Madori (|201 If ) for Spitzer 
wavelengths where the scatter about the TF relation is 
reduced from 0.43 in the extinction-corrected B— band 
data to 0.31 in the non-extinction-corrected 3.6/im data. 
One avenue for advancement is t herefore the use o f large 
IR photometric galaxy samples (|Sheth et al.ll2010T ). 

The most straightforward way for advancement is to 
obtain kinematic measurements for a large number of 
galaxies with existing distance measurements. In Fig- 
ure [7| I present a subset of galaxies with SNe la distance 
measurements and highlight those for which all the nec- 
essary data are available. An interesting indication of 
how these data may help refine distances is provided by 
two outliers from the FM (at 50 Mpc and the one da- 
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Figure 6. Hubble diagram using distances derived using the FM. 
Pressure supported systems are plotted as red circles and rota- 
tionally supported systems as blue squares. Inset shows data at 
v > 1500 km sec" 1 binned so that each bin contains 10 galaxies. 
The fitted line is constrained to go through the origin and corre- 
sponds to Ho = 78 ±2 (random) ±4 (systematic) km sec - 1 Mpc -1 
(adopted directly from Figure 12 of[Zaritsky ct al. 2 0111 ) 
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Figure 7. SN la distances and the FM. The upper panel shows 
the residuals about the FM relation (AFM) for galaxies with SN 
la distances. Blue and red symbols denote rotationally and pres- 
sure supported systems respectively. The bottom panel shows the 
Hubble diagram for galaxies with SN la distances available in the 
NED 1-D database. The galaxies with all the necessary data to 
be placed on the FM are highlighted in red. There are numer- 
ous galaxies with existing SN la distances that could be added to 
the FM with reasonable observational effort. The most noticeable 
outliers from the FM are also outliers in the Hubble diagram. 

turn beyond 100 Mpc), which are both outliers in the 
SNe Hubble diagram. That these galaxies are outliers 
in both panels suggests that there is a problem with the 
distance. With larger samples, one could then attempt 
to find the cause for the distance anomaly by comparing 
with other characteristics, either those of the SNe itself 
or those of the host galaxy. Refining the SNe distances 
has broad implications for cosmology. 



4.2. Using the Scaling Relationships for Mass 
Measurements 

The stellar mass of a galaxy is intricately con- 
nected to a variety of g alaxy properties: environ- 
ment (Kauffman n et ail I2004D . metal abundance 



Tremonti et al.l 20041) . star formation history 
Bundv et al. 20061) . and dark matter halo mass 
van den Bosch et^l2l)0l . ;o name a few. However, 
the stellar mass estimates we rely on are quite crude 
and potentially rife with systematic errors. A seemingly 
straightforward way to estimate the stellar mass is to 
use measurements of the stellar populations, such as 
color, to estimate T* and to convert the luminosity 
into a mass. However those estimates depend on 
three factors that are not well understood: 1) the 
galaxy's star formation history, which can be quite 



Harris fc Zarits kv 2004, 2009; 
Weisz et al.l 1201 ll) 2) stars' 



complex (see for example s 
lEskew fe Zaritikvl [20111: 

behavior during the phase(s) of their life at which 
they are at thei r most luminous, which is probl e matic 
dMaraston et al l [20061 iConrov. Gunn. fc White! [200l 
IMcQuinn. et al.l 120111: Melbourne et al.l 120121 ) and 3) 
the initial distribution of stellar masses (the initial mass 
function or IMF), which is a well-known unresolved 
problem (Bas tian et al.|[2010l) . An alternative method is 
to dynamically measure the mass and, if necessary, make 
a correction for the amount of dark matter present. The 
FM enables us to measure the total masses within r^. 

An interesting set of stellar systems to examine are 
stellar clusters, which presumably contain no dark mat- 
ter and are single age populations. As such, the es- 
timates of Y/i obtained using the FM provide a mea- 
surement of the mass-to-light ratio of a stellar pop- 
ulation, T„, of a certain age. A large compilation 
of Local Group stellar clust e r dat a was published by 
McLaughlin & van der Marcl (2005), although the veloc- 
ity dispersions available at the time of that study were 
either for the older (> 10 Gyr) clusters or the ve ry young 
ones (< 100 million yr). I Zaritskv etldl (|2012D filled in 
this range by observing clusters over the full range of 
ages. 

Discrepancies between models and observations of 
galaxies, when found, are often attributed to de- 
viations in the IMF from the adopt ed prescrip- 
tion ( for some recent example s see Ivan Dokku 
2008t IDabringhausen et all 120091: iTreu et aLJ 1201 



van Dokkum fc Conrovll2010l: iDutton et al.ll2012f) rather 



than problems with either the star formation history or 
stellar evolutionary models. However, direct measure- 
ments of the initial mass function are difficult for various 
reasons (see lBastian et ail l2010. for a review), particu- 
larly over the full range of environments and conditions. 
These clusters provide one of the most direct tests of the 
IMF. 

The princi pal empirical result presented by 

IZaritskv et al.l ()2012D . the relationship between T* 
and age, is reproduced in Figure |S] (their Figure 9 
reproduced here). The naive expectation, that T* 
rises continually with age is evidently not reproduced 
by the data. Internal dynamics, principally two-body 
relaxation which causes preferential loss of lo w mass 
stars over time (jSpitzerl Il958t IKruiissenl 120081) , plays 
a role here, but it is insufficient to explain the large 
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drop i n T* between a few and 10 Gyr. iZaritskv et al.l 
(2012) sup port this clai m by a pplying the dynamical 
models of lAnders et al.l (|2009[ ). A second possible 
reason for why the naive expectation may not be met 
is the influence of binary stars, which can artificially 
inflate the obser ved velocity dispersion of a system 
(for examples see iMinor etaJ 120101: iGieles et al.ll2010D . 
Again the effects are considered and found to be too 
weak to explain the observations or too contrived, 
requiring unknown an binary population that affects 
the kinematics of ~ 1 Gyr populations but not those of 
age > few Gyr. Variations in the initial mass function 
among these systems is a viable exp l anatio n. 

In the same Figure, IZaritskv et al.l (|2012l ) also included 
two model tracks for the evolution of simple stellar pop- 
ulations. First, the plotted values of T* for the stel- 
lar clusters are corrected for the effect o f internal dy- 
namic al evolution using the results of the lAnders et al.1 
2009) models. They plot the res ults of P EGASE models 



Fioc fc Rocca-Volmerangell997|) using a lSalpeterl (|1955[ ) 
IMF, spanning from 0.1 to 120 M©, with default stel- 
lar mass loss and binarity parameters, and metallicity 
matching the mean of the young clusters (—0.4). These 
model values are not renormalized in any way and yet 
do an acceptable job of reproducing the trend seen in 
the younger clusters, for clusters of 8 < log (age [yrs]) 
< 9.4. Discrepancies at younger ages can be ignored 
due to the likelihood that these clusters are not relaxed 
([Goodwin fe Bastian 2 00(f). T he second model is that of 
a lightweighted lKroupal (j200lD IMF, where ~ 50% of the 
mass is removed (lightweighted) to produce the match 
in the Figure. A iCharbrierl (|2003D IMF works similarly. 
One interpretation of this Figure hypothesizes the pres- 
ence of two IMFs for star clusters, one IMF being primar- 
ily, but not exclusively, appropriate for older, metal poor 
clusters and the other for primarily, but not exclusively, 
for younger, metal rich clusters. The young (log(age 
[yrs])<9.5) clusters are well-described by a bottom- heavy 
IMF, such as a Salpeter IMF, while the older clusters are 
better described by a top-heavy IMF, such as a light- 
weighted Kroupa IMF, although neither of these specific 
forms is a unique solution. Ongoing work will at least 
double the sample size, but we will eventually need to 
obtain high resolution deep imaging of the intermediate 
age clusters to confirm that these have a bottom heavy 
mass function. While some work has been done with 
HST (jGlatt et al.ll20lil ). deeper, less ambiguous results 
are necessary to definitively prove or disprove this inter- 
pretation. Such work is within reach of adaptive optics 
on large telescopes. 

In addition to stellar clusters, the FM can also be ap- 
plied to estimate the masses (within r/J of galaxies across 
cosmic time. In particular, this approach can test the hy- 
pothesis that large deviations from stellar evolutionary 
models, rather than IMF variations, are responsible for 
sudden drop in T» after a few Gyr. Using measurements 
of the ages and structural para meters of two indepen - 
dent sets of early type galaxies, IZaritskv et al.l (|2012p j 
compared the derived values of Th to those derived for 
the clusters. Such comparisons constrain both Y* and 
the fraction of the mass in the form of dark matter 
within r/j. They used both a study of Sloa n Digital Sky 
Survey (SDSS) local galaxies presented by iGraves et al.l 
(2009) and a study of galaxy cluster FP measurements 




Figure 8. The evolution of Y„ as inferred from early type galaxies 
and comparison to models and our cluster data. The stellar cluster 
data are plotted in solid circles (black), the results obt ained us- 
ing the results from van Dokkum & van dcr Marcl ( 2007) arc plot- 
ted with open circle s (red) when norma lized to match the results 
from the sample of [Graves et al. (20 Q3), which are plotted with 
square s (pink). Also shown are the Ivan Dokkum van der Marell 
(2007) data assuming a dark matter fraction of 80%, so as to qual- 
itatively match to the stellar cluster data, in open circles (blue). 
The dark matter fraction is likely to be somewhere between and 
80%. The upper solid line represents the values of T/ t from a 
PEGASE model of a population with an instantaneous burst at 
t = and a Salpeter IMF, while the lower represents a model with 
a light-weighted Kro upa IMF (adopted directly from Figure 9 of 
IZaritskv etaLl fWv3) ). 



presented by Ivan Dokkum fc van der Marell (f2007ri . The 
latter study provides only differences in Th as a 
function of redshift, so global shifts of are al- 
lowed. The results of those studies are include in Fig- 
ure [8] for comparison to the cluste r values (after the 
Ivan Dokkum fc van der Marell (120071) d ata were normal- 
ized to match the IGraves et alJ (|2009l ) data). An im- 
portant distinction for the galaxy results, is that these 
systems do contain dark matter and that the exact pro- 
portion of dark matter within is unknown and likely 
to vary as a function of the structural parameters (for 
some examples from the long history of this topic see 
i Babul fc Reed Il992t lGrahamll2002t iMarinoni fc Hudson! 
20021: Ivan den Bosch et al.1 12003 iCappellari et al.l 120061: 
IZaritskv et al.H2006aHWolf et al.ll2010D . 

The comparisons between clusters, galaxies and mod- 
els provide several interesting results. The early- type 
galaxies lie between the extrapolation of the young clus- 
ter trend and that of the old clusters. Because of the 
dark matter content of galaxies, the actual value of T* 
for the galaxies is likely to be lower than that plotted. 
This suggests that they will not lie on the track defined 
by the young clusters, which is already slightly above 
the galaxy datcQ- On the other hand, unless the dark 
matter fraction is quite large (~ 80%, the blue points in 
the Figure), the ellipticals will not resemble the old clus- 
ters either. Although the possibility of such a large DM 
fraction is not excluded, we may be seeing the effects of 
having a mixture of populations of stars, some with the 
IMF of the old clusters and some with that of the young 
clusters. However, there is no strong deviation from the 
general behavior predicted by simple stellar evolutionary 
models and so that does not appear to be the cause of 

6 Careful modeling is necessary to reach robust conclusions be- 
cause the galaxies are not necessarily single age populations. If an 
old galaxy has a small population of younger stars, those stars will 
lower the effective T of the galaxy. 
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the behavior of T* in the stellar cluster sample. 

5. SUMMARY 

Order implies rules. Rules governing the structure of 
dynamical systems are the manifestation of underlying 
physics. Understanding this physics is the overarching 
goal of the study of galaxy formation and evolution. Is 
there order among galaxies? 

In the early days of modern astronomy, the morpholog- 
ical appearance of galaxies was the principal characteris- 
tic that dominated att e mpts at descri bing order among 
galaxies (|Hubbld Il926t ISandagdli"96l . This approach 
has been mostly upended by 1) the realization that ac- 
cretion and mergers drive galaxy evolution and alter mor- 
phologies, 2) the dynamical importance of dark matter, 
which clearly is not a part of a morphological scheme, 
3) the increasingly quantitative nature of astronomy and 
large digital surveys, and 4) the discovery of scaling re- 
lations such as the Fundamental Plane (FP) and Tully- 
Fisher (TF). In this review, I have described the devel- 
opment of a scaling relation that is an extended version 
of the FP, called the Fundamental Manifold (FM), that 
works on all classes of stellar systems and exhibits a scat- 
ter comparable to the more restrictive FP and TF rela- 
tions. 

The FM relationship depends on three distinct condi- 
tions being satisfied: 1) that the virial theorem is ap- 
plicable, 2) that the derived mass enclosed within the 
half-light radius, r/,, be at most only weakly dependent 
on the distribution function of the tracer particles and 
the gravitational potential, and 3) that the mass-to-light 
ratio within rh, "T\, depend on at most the observed 
parameters Vh and 1^. The first is as expected be- 
cause we are working with systems that are virialized 
(r/, -C f20o)- The second is somewhat surprising, but 
has now been verified at least for spheroidal systems, 
which are the ones that are most likely to have large 
variations among d istri bution functions, b y the work of 
iWalker et all (|200l and I Wolf et all ([2010 ). The last still 
requires explanation, although it is manifestly a critical 
component of a complete theory of galaxy formation and 
evolution. 

Going forward there are several areas that will yield 
immediate returns. First, homogenous samples across 
galaxy type with detailed kinematics so that there is uni- 
form, high quality data across the parameter space are 
critical to testing Equation measuring the parameter 
that relates velocity dispersion and rotation, a, and de- 
termining the intrinsic scatter about the relation as a 
function of galaxy mass and morphology. Second, an im- 
proved understanding of the initial mass function is key 
to accurate and precise calculations of T*. Determin- 
ing T* is necessary to uncovering the behavior of dark 
matter as a function of mass and galaxy type. Third, 
direct application of the FM to a number of topics, in- 
cluding distance determinations, will yield results even 
if a deeper understanding of galaxy formation remains 
elusive in the short term. The FM, like the FP and TF 
before it, should help propel research forward in a wide 
range of topic areas. 

The author acknowledges the close collaboration with 
Ann Zabludoff and Anthony Gonzalez during which most 



of the material reviewed here was developed. The au- 
thor acknowledges financial support for this work from a 
NASA LTSA award NNG05GE82G and NSF grant AST- 
0307482 as well as the hospitality of NYU CCPP during 
his many visits. 

REFERENCES 



Aaronson, M., & Mould, J. 1983, ApJ, 265, 1 
Anders, P., Lamers, H.J.G.L.M., & Baumgardt, H. 2009, A&A, 
502, 817 

Anderson, M.E., & Bregman, J.N., 2010, ApJ, 714, 320 
Babul, A., & Rees, M.J. 1992, MNRAS, 255, 346 
Bahcall, J. N., & Tremaine, S. 1981, ApJ, 244, 805 
Bastian, N., Covey, K.R., & Meyer, M.R. 2010, Ann. Rev. 

Astron. Astroph., 48, 339 
Bedregal, A.G., Aragon-Salamanca, A., & Merrifield, M.R. 2006, 

MNRAS373, 1 125 
Bender, R. 1990, A&A, 229, 441 
Bernardi, M. et al. 2003, AJ, 125, 1866 
Bernardi, M., Roche, N, Shankar, F., & Sheth, R.K. 2011, 

MNRAS, 412, 6 
Bolton, A. S., Treu, T., Koopmans, L. V. E., Gavazzi, R., 

Moustakas, L. A., Buries, S., Schlegel, D. J. & Wayth, R., 2008, 

ApJ, 683, 248 
Bundy, K. et al. 2006, ApJ, 651, 120 

Burstein, D. & Bender, R. & Faber, S. & Nolthenius, R., 1997, 

AJ, 114, 1365 
Cappellari, M. et al. 2006, MNRAS, 366,1126 

Cappcllari, M. et al. 2012, submitted to MNRAS. larXiv:120 8. 3523 

Chabrier, C, 2003, PASP, 115, 763 

Chiboucas, K. et. al. 2011, ApJ, 737, 86 

Conroy, C, Gunn, J. E. & White, M., 2008, 

Courteau, S., Dutton, A. A., van den Bosch, F.C., MacArthur, 

L.A., Dekel, A., Mcintosh, D.H., & Dale, D.A. 2007, ApJ, 671, 

203 

Cretton, N, de Zeeuw, P.T., van der Marel, R.P., & Rix, H.W.R. 

1999, ApJ, 124, 383 
Dabringhausen, J., Kroupa, P., Baumgardt, H. 2009, MNRAS, 

394, 1529 

Dalcanton, J. J., Spergel, D.N., & Summers, F.J. 1997, ApJ, 482, 
659 

Davies, R. L., Efstathiou, C, Fall, S.M., Illingworth, C, 

Schechter, P.L. 1983, ApJ, 266, 41 
Dejonghe, H. 1987, MNRAS, 133, 217 
Djorgovski, S. & Davis, M., 1987, ApJ, 313, 59 
Dressier, A. & Lynden-Bell, D. & Burstein, D. & Davies, R. L., 

Faber, S. M., Terlevich, R. & Wegner, C, 1987, ApJ, 313, 42 
Drinkwater, M.J., Jones, J.B., Gregg, M.D., & Phillipps, S. 2000, 

PASA, 17, 227 
Dutton, A.A., 2012, MNRASin press (arXiv: 1206: 1855) 
Dutton, A.A., Mendel, J.T., & Simard, L. 2010, MNRAS, 422, 33 
Dutton, A. A. & van den Bosch, F.C. 2012, MNRAS, 421, 608 
Erickson, K., Gottesman, S.T., & Hunter, J.H. Jr., 1987, Nature, 

325, 779 

Eskew, M. & Zaritsky, D. 2011, AJ, 141, 69 

Eskew, M., Zaritsky, D., & Meidt, S. 2012, AJ, 143, 139 

Faber, S.M. & Jackson, R.E. 1976, ApJ, 204, 668 

Fall, S. M., & Efstathiou, G. 1980 MNRAS, 193, 189 

Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 

Forbes, D. A., Lasky, P., Graham, A. W., & Spitler, L., 2008, 

MNRAS, 389, 1924 
Freedman, W.L., & Madore, B.F., 2011, ApJ, 734, 46 
Gallagher, J. S., Garnavich, P. M., Caldwell, N., Kirshner, R. P., 

Jha, S. W., Li, W., Ganeshalingam, M. and Filippenko, A. V., 

2008, ApJ, 685, 752 
Gieles, M., Sana, H., & Pertegies Zwart, S.F. 2010, MNRAS, 402, 

1750 

Glatt, K. et al. 2011, AJ, 142, 36 

Goodwin, S. P., & Bastian, N. 2006, MNRAS, 373, 752 

Graves, G.J., Faber, S.M., & Schiavon, R.P. 2009, ApJ, 698, 1590 

Graham, A. 2002, MNRAS, 334, 721 

Harris, J., & Zaritsky, D. 2004, AJ, 127, 1531 

Harris, J., & Zaritsky, D. 2009, AJ, 138, 1243 

Heisler, J., Tremaine, S., & Bahcall, J.N. 1985, ApJ, 298, 8 



12 



Zaritsky 



Hilker, M., Infante, L., Vieira, G., Kissler-Patig, M., & Richtler, 

T. 1999, A&AS, 134, 75 
Hoekstra, H., Hsieh, B.C., Yee, H.K.C., Lin, H., & Gladders, 

M.D. 2005, ApJ, 635, 73 
Howell, D.A. et al. 2009, ApJ, 691, 661 
Hubble, E.P. 1926, ApJ, 64, 32 

J0rgensen, I., Franx, M., & Kjaergaard, P. 1996, MNRAS, 280, 167 
Kassin, S.A., et al. 2007, 660, 35 

Kaufmann, G., White, S. D. M., Heckman, T., Menard, B, 

Brinchmann, J., Chariot, S., Tremonti, C, & Brinkmann, J. 

2004, MNRAS, 353, 713 
Khosroshahi, H.G., Wadadekar, Y., Kcmbhavi, A. & Mobasher, 

B. 2000, ApJ, 531, 103 
Kormendy, J. 1985, ApJ, 295, 73 
Kroupa, P. 2001, MNRAS, 322, 231 
Kruijssen, J. M.D. 2008, A&A, 486, 21 
Little, B., & Tremaine, S. ApJ, 320, 493 
Mandelbaum, R., Seljak, U., Kauffmann, G., Hirata, C, & 

Brinkmann, J. 2006, MNRAS, 368, 715 
Maraston, C, Daddi, E., Renzini, A, Cimatti, A. & Dickinson, 

M., Papovich, C, Pasquali, A., & Pirzkal, N., 2006, ApJ, 652, 

85 

Marinoni, C, & Hudson, M.J. 2002, ApJ, 569, 101 
McGaugh, S.S, 2012 AJ, 143, 40 

McLaughlin, D. E. & van der Marel, R. P., 2005, ApJS, 161, 304 
McQuinn, K.B.W., Skillman, E.D., Dalcanton, J. J., Dolphin, 

A.E., Holtzman, J., Weisz, D., & Williams, B.F. 2011, AJ, 740, 

48 

Mcidt, S. et al. 2012, ApJ, 744, 17 
Melbourne, J., et al. 2012, ApJ, 748, 47 
Merrifield, M.,R. & Kent, S.M. 1990, ApJ, 99, 1548 
Mieske, S. et al. 2008, A&A, 487, 921 

Minor, Q.E., Martinez, G., Bullock, J., Kaplinghat, M., & 

Trainor, R. 2010, ApJ, 721, 1142 
Misgeld, I. & Hilker, M. 2011, MNRAS, 414, 3699 
Misgeld, I., Mieske, S., Hilker, M., Richtler, T., Gerogiev, I.Y. & 

Schuberth, Y. 2011, A&A, 531, 4 
Mould, J., Mingsheng, H., & Bothun, G. 1989, ApJ, 347, 112 
Page, T. 1952, ApJ, 116, 63 
Pasquiato, M. & Bertin, G. 2008, a, 489, 1079 
Persic, M., & Salucci, P. 1991, MNRAS, 248, 325 
Phillipps, S., Drinkwater, M.J., Gregg, M.D., & Jones, J.B. 2001, 

ApJ, 560, 201 
Pizagno, J., et al. 2007, AJ, 134, 945 
Romaniello, M. et al. 2008, A&A, 488, 731 
Saintonge, A., & Spckkcns, K. 2011, ApJ, 726, 77 
Salpeter, E.E. 1955, ApJ, 121, 161 

Salucci, P., Lapi, A., Tonini, C, Gentile, G., Yegorova, I., & 
Klein, U. 2007, MNRAS, 378, 41 



Salucci, P., Wilkinson, M.I., Walker, M.G., Gilmore, G.F., 

Grebel, E.K., Koch, A., Frigerio Martins, C, & Wyse, R. F. G. 
2012, MNRAS, 420, 2034 

Sandage, A. 1961, The Hubble Atlas of Galaxies, (Washington: 
Carnegie Institution) 

Shcth, K., et al. 2010, PASP, 122, 1397 

Smith, M.A. & Gray, D.F. 1976, PASP, 88, 809 

Spitzer, L. 1958, ApJ, 127, 17 

Springel, V. et al. 2005 Nature, 435, 629 

Statler, T.S. 1987, ApJ, 321, 113 

Tollerud, E.J., Boylan-Kolchin, M., Barton, E.J., Bullock, J.S., & 

Trihn, C.Q. 2011, ApJ, 738, 102 
Tollerud, E.J., Bullock, J.S., Strigari, L. E., & Willman, B. 2008, 

ApJ, 688, 277 
Tremonti, C.A. et al. 2004, ApJ, 613, 898 
Treu, T., Auger, M.W., Koopmans, L.V.E., Gavazzi, R., 

Marshall, P.J., & Bolton, A.S. 2010, ApJ, 709, 1195 
Tully, R. B. and Fisher, J. R., 1977, A&A, 54, 661 
van den Bosch, F. C, Yang, X., Mo, H.J. 2003, MNRAS, 340,771 
van der Marel, R.P., & van Dokkum, P.G. 2007, ApJ, 668, 756 
Vandervroot, P. O. 1984, ApJ, 287, 475 
van Dokkum, P.G., 2008, ApJ, 674, 29 
van Dokkum, P.G., & Conroy, C. 2010, Nature, 468, 940 
van Dokkum, P. G., & van der Marel, R. P. 2007, ApJ, 655, 30 
Walker, M.G., Mateo, M., Olszewski, E.W., Penarrubia, J., Wyn 

Evans, N., & Gilmore, G. ApJ, 704, 1274 
Weiner, B. et al. 2006, ApJ, 653, 1049 
Weisz, D., et al. 2011, ApJ, 739, 5 
Willman, B. & Strader, J. 2012, AJ, in press 

Wolf, J., Martinez, G.D., Bullock, J.S., Kaplinghat, M., Geha, M., 

Munoz, R.R., Simon, J.D., & Avedo, F.F. MNRAS, 406, 1220 
Yang, X., Mo, H.J., Jing, Y.P., & van den Bosch, F.C. 2005, 

MNRAS, 358, 217 
Yegorova, LA. & Salucci, P. 2007, MNRAS, 377, 507 
Zaritsky, D., Colucci, J., Pessev, P., Bernstein, R.A., & Chandar, 

R. 2012, ApJ, in press 
Zaritsky, D., Gonzalez, A. H. & Zabludoff, A. I., 2006a, ApJ, 638, 

725 

Zaritsky, D., Gonzalez, A. H. & Zabludoff, A. I., 2006b, ApJ, 642, 
L37 

Zaritsky, D. & White, S.D.M. 1994, ApJ, 435, 599 
Zaritsky, D., Zabludoff, A. I. & Gonzalez, A. H., 2008, ApJ, 682, 
68 

Zaritsky, D., Zabludoff, A. I. & Gonzalez, A. H., 2011, ApJ, 727, 
116 

Zaritsky, D., Zabludoff, A. I. & Gonzalez, A. H., 2012, ApJ, 758, 
15 

Zhu, G., Blanton, M.R., & Moustakas, J. 2010, ApJ, 722, 491 



